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Abstract 

The production of photons at very small angles with respect to the proton beam direction 
is studied in deep-inelastic positron-proton scattering at HERA. The data are taken with 
the HI detector in the years 2006 and 2007 and correspond to an integrated luminosity of 
126 pb^^. The analysis covers the range of negative four momentum transfer squared at 
the positron vertex 6 < < 100 GeV^ and inelasticity 0.05 < y < 0.6. Cross sections 
are measured for the most energetic photon with pseudorapidity rj > 7.9 as a function 
of its transverse momentum p!^"''^ and longitudinal momentum fraction of the incoming 
proton x^""^. In addition, the cross sections are studied as a function of the sum of the 
longitudinal momentum fraction xf^"^ of all photons in the pseudorapidity range r] > 7.9. 
The cross sections are normalised to the inclusive deep-inelastic scattering cross section 
and compared to the predictions of models of deep-inelastic scattering and models of the 
hadronic interactions of high energy cosmic rays. 
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1 Introduction 



Measurements of particle production at very small angles with respect to the proton beam di- 
rection (forward direction) in positron-proton collisions are important for the understanding 
of proton fragmentation. These measurements also provide important constraints for the mod- 
elling of the high energy air showers and thereby are very valuable for the understanding of high 
energy cosmic ray data. The HI and ZEUS experiments at the e^p collider HERA have pub- 
lished several analyses on the production of forward protons and neutrons which carry a large 
fraction of the longitudinal momentum of the incoming proton [UJ-15|. These measurements 
probe different mechanisms related to the baryon production in forward direction, such as elas- 
tic scattering of the proton, diffractive dissociation, pion exchange and string fragmentation. In 
particular, these measurements test the hypothesis of limiting fragmentation |l6ll3, according 
to which, in the high-energy limit, the cross section for the inclusive production of particles in 
the target fragmentation region is independent of the incident projectile energy. This hypothesis 
implies, that in deep-inelastic scattering (DIS) forward particle production cross sections are 
independent of the Bjorken-x and the virtuality of the exchanged photon Q"^. 

The measurement of the photon production in the forward direction can provide new input 
to the understanding of proton fragmentation, and is complementary to forward baryon mea- 
surements. The production of photons and tt^ mesons in the proton fragmentation region has 
been studied in pp and pp collisions at SPS and the LHC colliders (Hllll. The analysis pre- 
sented here is the first measurement of very forward photons in DIS e+p collisions at HERA. 
The photons are detected at very small angles below 0.75 mrad with respect to proton beam 
direction. It relies on the upgraded HI Forward Neutron Calorimeter (FNC) which includes an 
electromagnetic section. 



2 Experimental Procedure and Data Analysis 

The data used in this analysis were collected with the HI detector at HERA in the years 2006 and 
2007 and correspond to an integrated luminosity of 126 pb^^. During the period corresponding 
to the analysis data set HERA collided positrons and protons with energies of E^, = 27.6 GeV 
and Ep = 920 GeV, respectively, corresponding to a centre-of-mass energy of ^/s = 319 GeV. 

2.1 HI detector 

A detailed description of the HI detector can be found elsewhere fTOj-fTSl. Only the detector 
components relevant to this analysis are briefly described here. The origin of the right-handed 
HI coordinate system is the nominal e~^p interaction point. The direction of the proton beam 
defines the positive z axis; the polar angle 9 is measured with respect to this axis. Transverse 
momenta are measured in the x-y plane. The pseudorapidity is defined by = — In (tan |) 
and is measured in the lab frame. The polar angles 9 < 0.75 mrad correspond to pseudorapidity 
range rj > 7.9. 
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The interaction region is surrounded by a two-layer silicon strip detector and two large 
concentric drift chambers. Charged particle momenta are measured in the angular range 
25° < 9 < 155°. The tracking system is surrounded by a finely segmented Liquid Argon (LAr) 
calorimeter, which covers the polar angle range of 4° < 6* < 154° with full azimuthal accep- 
tance. The LAr calorimeter consists of an electromagnetic section with lead absorber and a 
hadronic section with steel absorber. The total depth of the LAr calorimeter ranges from 4.5 
to 8 hadronic interaction lengths. The backward region (153° < 9 < 177.8°) is covered by a 
lead/scintillating-fibre calorimeter (SpaCal). Its main purpose is the detection of the scattered 
positron. The energy resolution for positrons is a{E)/E ^ 7.1%/ 



^^[GeV] © 1%, as deter- 
mined in test beam measurements [[T4ll . The LAr and SpaCal calorimeters are surrounded by 
a superconducting solenoid which provides a uniform magnetic field of 1.16 T along the beam 
direction. 

The luminosity is measured via the Bethe-Heitler Bremsstrahlung process ep — j- e'p'j, the 
final state photon being detected in a tungsten/quartz-fibre sampling calorimeter at z = — 103 m. 

The data sample of this analysis was collected using triggers which require the scattered 
positron to be measured in the SpaCal. The trigger efficiency is about 96% for the analysis 
phase space as determined from data using independently triggered data. 



FNC 

MAIN CALORIMETER 




(b) 



PRESHOWER 




Figure 1: (a) A schematic view of the HI FNC. (b) Layout of 9 vertical and 9 horizontal read- 
out strips of the Preshower Calorimeter. The hatched area shows the geometrical acceptance 
window defined by the beam-line elements. The area corresponding to r] > 7.9 is indicated by 
dashed circle. 



2.2 Detection of forward neutral particles 

Neutral particles produced at very small polar angles can be detected in the FNC calorimeter, 
which is situated at a polar angle of 0°, at 2; = +106 m from the interaction point. A schematic 
view of the HI FNC used during the HERA-II running period is shown in figure [T^. A detailed 
description of the detector is given in [5]. The FNC consists of the Main Calorimeter and the 
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Preshower Calorimeter. The Main Calorimeter is a lead-scintillator sandwich calorimeter with 
a total length of 8.9 nuclear interaction lengths. The Preshower Calorimeter is a 40 cm long 
lead-scintillator sandwich calorimeter. The length corresponds to about 60 radiation lengths. 
The Preshower Calorimeter is composed of 24 planes: the first 12 planes each consist of a lead 
plate of 7.5 mm thickness and a scintillator plate of 2.6 mm thickness. The second 12 planes 
each consist of a lead plate of 14 mm thickness and a scintillator plate of 5.2 mm thickness. 
The transverse size of the scintillating plates is 26 x 26 cm^. Each scintillating plate has 45 
grooves with 1.2 mm wavelength shifter fibres attached down one side. The orientation of 
fibres alternates from horizontal to vertical in consecutive planes. For each plane, the fibres are 
bundled into nine strips of five fibres each. Longitudinally, all strips are combined leading to 9 
vertical and 9 horizontal towers which are finally connected to 18 photomultipliers. 

The acceptance of the FNC is defined by the aperture of the HERA beam-line magnets and 
is limited to scattering angles of 6* £ 0.8 mrad with approximately 30% azimuthal coverage. The 
geometrical acceptance window of the FNC is shown in figure [IJ5 together with the layout of 
the Preshower Calorimeter readout strips. 

The longitudinal segmentation of the FNC allows efficient discrimination of photons from 
hadrons. The photon reconstruction algorithm is based on the fact that electromagnetic showers 
are fully contained in the Preshower Calorimeter with no energy deposits above the noise level 
in the Main Calorimeter. For high energy neutrons most of the energy is contained in the 
Main Calorimeter. However, low energy neutrons deposit large fractions of their energy in the 
Preshower Calorimeter. The fraction of neutrons which can be misidentified as photons is about 
10% for 90 GeV neutrons decreasing to below 1% for neutrons with an energy of 200 GeV, 
as determined from the Monte Carlo (MC) simulation. The energy deposits in the FNC which 
are contained in the Preshower Calorimeter are classified as electromagnetic clusters and are 
considered as photon candidates. The detection and reconstruction efficiency for photons in the 
measured angular range 9 < 0.75 mrad, as estimated from MC simulation, is about 85% for 
100 GeV photons increasing to 95% for photons with energies of 900 GeV. Losses are mainly 
due to interactions with the beampipe. 

All modules of the FNC were initially calibrated at CERN using 120-230 GeV electron 
and 120-350 GeV hadron beams. After the calorimeter was installed at DESY, the stability of 
calibration constants was monitored using interactions between the proton beam and residual 
gas in the beam pipe, as described in Refined calibration constants for electromagnetic 
showers are determined using an iterative procedure based on the assumption that the maximum 
photon energy, EJ^"'^, as measured in the Preshower Calorimeter, is expected to be equal to the 
proton beam energy in case of unlimited statistics. This calibration procedure also utilises data 
from HERA runs with reduced proton beam energies of 460 GeV and 575 GeV. The validity of 
this algorithm is tested with MC simulation. 

The measured photon energy spectra for the three proton beam energies are displayed in 
figure [2^. The correlation between the beam energy and the maximum photon energy i^^"^ as 
determined by the iterative procedure and after applying the calibration is shown in figure [Ab- 
using this calibration procedure, the linearity of the energy response and the absolute energy 
scale are verified to a precision of 5%. 
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EFNc[GeV] Ep[GeV] 

Figure 2: (a) The measured photon energy spectra for three proton beam energies, (b) The 
correlation between the proton beam energy and E!^"-^. 



The energy resolution of the FNC calorimeter for electromagnetic showers is a{E)/E ^ 
20%/ ^/E~\^eV] © 2%, as determined in test beam measurements. The spatial resolution for 
single electromagnetic showers and for those hadronic showers which started to develop in the 
Preshower Calorimeter is about 2 mm. 



2.3 Kinematics and event selection 

The kinematic variables used to describe high energy DIS interactions are the exchanged photon 
virtuality Q^, the inelasticity y and the Bjorken scaling variable xsj- They are defined as 

Q = -g , XBj = , y = r , (1) 

2p ■ q p ■ k 

where p, k and q are the four-momenta of the incident proton, the incident positron and the 
virtual photon, respectively. 

The selection of DIS events is based on the identification of the scattered positron as the most 
energetic compact calorimetric deposit in the SpaCal with an energy E'^ > 11 GeV and a polar 
angle 156° < 9'^ < 175°. The z coordinate of the primary event vertex is required to be within 
±35 cm of the nominal position of the interaction point. The hadronic final state is reconstructed 
using an energy flow algorithm which combines charged particles measured in the tracker with 
information from the SpaCal and LAr calorimeters Ill7[|18i] . To suppress events with hard initial 
state radiation, as well as events originating from non-e+p interactions, the quantity E — p^, 
summed over all reconstructed final state particles including the positron, is required to lie 
between 35 GeV and 70 GeV. This quantity, which uses the energy and longitudinal momentum 
component of each final state particle, is expected to be twice the electron beam energy for fully 
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contained events. In addition, events are restricted to the kinematic range 6 < < 100 GeV^ 
and 0.05 < y < 0.6. These variables are reconstructed using a technique which optimises 
the resolution throughout the measured y range, exploiting information from both the scattered 
positron and the hadronic final state [16]. The data sample of the DIS events contains about 
9.2 million events. 

Events containing forward photons are selected by requiring an electromagnetic cluster in 
the FNC with a pseudorapidity above 7.9 and an energy above 92 GeV, which corresponds to 
the longitudinal momentum fraction xl = E^/Ep > 0.1, where Ep and E^y are the proton beam 
and forward photon energy, respectively. The data sample contains 78740 events. 

In this analysis normalised differential cross sections are measured for the most energetic 
forward photon (leading photon) with pseudorapidity 7] > 7.9 as a function of its longitudinal 
momentum fraction x^"'' and transverse momentum p?^'"^, in the range 0.1 < < 0.7. Cross 
sections are also measured as a function of the sum of longitudinal momentum fractions of all 
forward photons with t] > 7.9, xf"™ = in the range 0.1 < xf"™ < 0.95. These cross 

sections are given as the fraction of DIS events having forward photon in the r]-XL regions 
given above. Finally, the ratio of the forward photon production cross section to the inclusive 
DIS cross section is presented as a function of Q"^ and xbj- 

2.4 Monte Carlo simulations and corrections to the data 

Monte Carlo simulations are used to correct the data for the effects of detector acceptance, in- 
efficiencies, migrations between measurement bins due to finite detector resolution and QED 
radiation from the positron. All generated events are passed through a GEANT3 [fT9ll based sim- 
ulation of the HI apparatus and are then processed using the same reconstruction and analysis 
chain as is used for the data. 

The DJANGOH 1203 program is used to generates inclusive DIS events. It is based on 
leading order electroweak cross sections and takes into account QCD effects up to order a^. 
Higher order QCD effects are simulated using leading log parton showers as implemented in 
LEPTO [21], or using the Colour Dipole Model (CDM) as implemented in ARIADNE [i22i. 
Subsequent hadronisation effects are modelled using the Lund string fragmentation model as 
implemented in JETSET Il23ll24ll . Higher order electroweak processes are simulated using an 
interface to HERACLES [|25l . The LEPTO program includes the simulation of soft colour 
interactions (SCI) [|26ll . in which the production of diffraction-like configurations is enhanced 
via non-perturbative colour rearrangements between the outgoing partons. In the measured xl 
range, omitting the SCI in the LEPTO would decrease the predicted yield of forward photons 
by 5% at lowest xl and 2% at highest xl- The simulation with CDM uses the parameters tuned 
to describe HI forward jet measurements [[27|. The DJANGOH MC simulations are calculated 
using the HIPDF 2009 parameterisation [28] of the parton distributions in the proton. In the 
following, the predictions based on LEPTO and ARIADNE are denoted LEPTO and CDM, 
respectively. 

The measurements are also compared with the predictions of several hadronic interac- 
tion models which are commonly used for the simulation of cosmic ray air shower cascades: 
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EPOS [|29l, QGSJET 01 Il30l[311. QGSJET II ||32ll33l and SIBYLL [|34l[35l. These phenomeno- 
logical models, based on general principles such as unitarity and analyticity of scattering am- 
plitudes, are combined with perturbative QCD predictions for high-p^- processes to obtain a 
description of the final state. The programs are interfaced with the PHOJET program [,36,1 for 
the simulation of e+p interactions. 

In all of these models, the main source of forward photons is the decay of 7r° mesons pro- 
duced from the hadronisation of the proton remnant. The measured distributions may contain 
background arising from several sources. The background from photoproduction processes, 
where the positron is scattered into the backward beam-pipe and a particle from the hadronic 
final state fakes the positron signature in the SpaCal, is estimated using the PHOJET MC gen- 
erator and found to be negligible. The selected sample may contain background from neu- 
trons reconstructed as electromagnetic clusters as explained above. For cluster energies above 
92 GeV this background is found to be negligible according to the MC simulation. The back- 
ground from the random coincidences of DIS events with a beam-related background signals 
in the FNC is estimated by combining DIS events with forward particles in adjacent bunch- 
crossings. It is found to be smaller than 1%. The background contributions are not subtracted 
from themeasured cross sections. 

Two or more particles entering the FNC are reconstructed as a single cluster due to the 
relatively large size of the FNC readout modules in combination with a small geometrical ac- 
ceptance window. According to the MC simulation, low energetic clusters reconstructed in the 
FNC mainly originate from single photons. The contribution from two photons increases almost 
linearly from 10% at about 450 GeV to 80% at 900 GeV (the contribution from three and more 
photons is below 1%). Therefore, the measurement of the cross section of single photon pro- 
duction is limited to xl < 0.7, while the measurement of the total forward photon production 
cross section is extended to larger xl- 

Factors determined from MC are used to correct distributions at the level of reconstructed 
particles back to the hadron level on a bin-by-bin basis. These correction factors include the 
effects of QED radiation from the positron. For the calculation of the correction factors the 
simulations are reweighted to describe the xl distributions of the data. The average of the 
correction factors determined from LEPTO and CDM is used. The size of the correction factors 
varies between 2 and 3.5 for x''l"'^, between 3 and 4 for x^"™, between 2.5 and 12 for p^"'^ 
and are about 3.2 for the and xbj distributions. They are dominated by the non-uniform 
azimuthal acceptance of the FNC, which is about 30% on average. The bin purities, defined 
as the fraction of events reconstructed in a particular bin that originate from that bin on hadron 
level, vary between 75% and 95%. 

2.5 Systematic uncertainties 

The systematic uncertainties on the cross section measurements are determined using MC sim- 
ulations, by propagating the corresponding uncertainty sources through the full analysis chain. 

As the cross sections are normalised to the inclusive DIS cross section measured in this ana- 
lysis, some important systematic uncertainties, such as the trigger efficiency, the luminosity and 
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the uncertainties related to the reconstruction of the scattered positron and of the hadronic final 
state are largely reduced or cancel. Uncertainties on the measurements of the scattered positron 
energy (1%) and angle (1 mrad), the energy of the hadronic final state (4%), and the uncertainty 
on the trigger efficiency (1%) lead to an average combined uncertainty of up to 2%. 

The absolute electromagnetic energy scale of the FNC is known to a precision of 5% as 
described in section [Z2l This leads to an uncertainty of 1% on the cross section measurement at 
low energies, increasing to 35% for the largest xl values. The acceptance of the FNC calorime- 
ter is defined by the interaction point and the geometry of the HERA magnets and is determined 
using MC simulations. The uncertainty of the impact position of the photon on the FNC is due 
to beam inclination and the uncertainty on the FNC position. It is estimated to be 5 mm. This re- 
sults in uncertainties on the FNC acceptance determination of up to 15% for the distributions 
and up to 60% for the p'^"'^ distribution. These effects are strongly correlated between measure- 
ment bins. For the and xbj measurements, these effects lead to normalisation uncertainty of 
approximately 7%. 

The systematic uncertainty arising from the model dependence of the data correction is taken 
as the difference of the corrections calculated using the LEPTO and CDM models. The resulting 
uncertainty on the cross-section increases from 1% to 6% for the x^"'' and pi^"''^ distributions, 
from 2% to 20% for the distribution, and from 1% to 2% for the and xbj distributions. 
Using different parton distribution functions in the MC simulation results in a negligible change 
in the cross section. 

The systematic errors shown in the figures and table are calculated as the quadratic sum 
of all contributions, which may vary from point to point. The total systematic error for the 
normalised cross section measurements ranges between 8% and 18% for x^"'^, 6% and 58% for 
p^"'^, 8% and 44% for ™ and 7% and 8% for and xbj- 

3 Results 

The measured normalised differential cross sections for the production of very forward photons 
in the pseudorapidity range r] > 7.9 in DIS in the kinematic range 6 < < 100 GeV^ and 
0.05 < y < 0.6, are presented in table [Hand figures [3ll5l The measurements are presented 
in figures [3] and |4] as a function of x^"'^ and p^^'"^ of the most energetic photon with 0.1 < 
^lead ^ Q 7 Theresuks as a function of the sum of longitudinal momentum fractions x|"™ of 
all photons with r] > 7.9 are presented in figure [51 

The data are compared with the predictions of models for inclusive DIS (LEPTO and CDM) 
and models of hadronic interactions (EPOS, SIBYLL and two versions of QGSJET). The ratios 
of MC model predictions to the measurements are shown separately. 

All models tested in this paper overestimate the total rate of forward photons. The LEPTO 
and CDM models predict about 70% more photons than measured, while EPOS, SIBYLL and 
QGSJET overestimate the rate of photons by about 30% to 50%. In contrast to the excess of 
photons in the CDM, the same model predicts a too low rate of forward neutrons as observed in 
previous HI analysis of forward neutron production [[Si- 
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The shapes of all measured distributions are well described by LEPTO. The CDM predicts 
harder xl and px spectra. The QGSJET model predicts slightly softer spectra in and px- The 
EPOS and SIBYLL models predict harder xl spectra, but describe reasonably the shape of px 
distribution. 

A measurement of the energy spectra of single photons produced at > 8.8 in pp collisions 
at 7 TeV centre-of-mass energy at the LHC has been recently reported by the LHCf Collab- 
oration BH. This measurement also shows significant discrepancies with the predictions from 
hadronic interaction models. However, a direct comparison of the HI and LHCf results is not 
possible due to the different kinematic ranges of the two measurements. 

The measurement of forward photons allows a test of the limiting fragmentation hypothesis, 
according to which the production of forward photons in DIS is insensitive to and xsj- 
To investigate this prediction, the ratio of the forward photon production cross section to the 
inclusive DIS cross section is measured as a function of and XBj (table |2] and figure [6l). 
Within the uncertainties the fraction of DIS events with forward photons is independent from 

and xsj in agreement with the limiting fragmentation hypothesis. A similar conclusion was 
obtained in the earlier HI analysis of forward neutron production [5J. The LEPTO and CDM 
predictions also included in figure [6] display a significant difference in normalisation compared 
to data as well as a slight dependence as a function of and XBj- 

4 Summary 

The production of high energy forward photons in the pseudorapidity range r] > 7.9 is studied 
for the first time at HERA in deep-inelastic positron-proton scattering in the kinematic region 
6 < < 100 GeV^, 0.05 < y < 0.6. The normalised DIS cross sections are presented for 
the production of the most energetic photon as a function of the longitudinal momentum frac- 
tion and transverse momentum in the range 0.1 < x^"'^ < 0.7, and as a function of the sum 
of longitudinal momentum fractions of all forward photons in the range 0.1 < xf"*" < 0.95. 
Predictions of Monte Carlo models overestimate the rate of photons. The shapes of the mea- 
sured cross sections are well described by the LEPTO MC simulation, while the colour dipole 
model predicts harder spectra in x^^ and pt- The measurement is also compared to predictions 
of models which are commonly used for the simulation of cosmic ray air shower cascades. All 
these models predict different spectra in xl and pt- None of the models can describe data in 
rate and in shape. Within the measured kinematic range, the relative rate of forward photons in 
DIS events is observed to be independent of and xbj, in agreement with the hypothesis of 
limiting fragmentation. The present measurement provides new information to further improve 
the understanding of proton fragmentation in collider and cosmic ray experiments. 
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ted sys. unc 
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0.00027 
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plead j-^j^gg 

[GeV] 


1 da 
^ DIS dpji^^ 
[GeV-i] 


^stat. 

[GeV-1] 


^ total sys. 

[GeV-1] 


^uncorrel.sys. 

[GeV-1] 


correla 

^Efnc 

[GeV-1] 


ted svs unc 

^XYfnc 
[GeV-1] 


3rtaintv 

^model 

[GeV-1] 


0.0^0.1 
0.1 -=-0.2 
0.2 H- 0.3 
0.3 -=-0.4 


0.159 
0.0971 
0.0220 
0.00395 


0.001 
0.0010 
0.0005 
0.00029 


0.010 
0.0116 
0.0056 
0.00229 


0.003 
0.0041 
0.0010 
0.00025 


0.005 
0.0068 
0.0024 
0.00087 
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0.0078 
0.0048 
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0.000468 
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0.00021 
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0.000039 
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0.000140 
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0.00041 

0.000070 


0.002 
0.0007 
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0.00031 

0.000050 



Table 1 : The normalised cross sections for the production of forward photons in the pseudora- 
pidity range 77 > 7.9 in deep-inelastic scattering in the kinematic region 6 < < 100 GeV^ 
and 0.05 < y < 0.6 as a function of the longitudinal momentum fraction x^'"' and transverse 
momentum p?^"'^ of the most energetic photon in the energy range 0.1 < x^"'^ < 0.7 and as a 
function of the sum of the longitudinal momentum fractions of photons a;£"™. For each mea- 
surement, the statistical, the total systematic, the uncorrected systematic uncertainties, and the 
bin-to-bin correlated systematic uncertainties due to the FNC absolute energy scale, the impact 
position of the FNC and the model dependence of data correction are given. 
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Q2 range [GeF^] 


^Dism 


^stat. 


^ total sys. 


Suncorrel.sys. 


correlated sys. uncertainty 






^model 


6.0- 
24.8- 
43.6- 
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- 62.4 
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0.0276 
0.0265 
0.0265 
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0.0003 
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0.0011 


0.0020 
0.0020 
0.0020 
0.0020 
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0.0003 
0.0003 
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0.0005 
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0.0011 
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0.0276 


0.0007 
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0.0006 


0.0011 


0.0017 
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Table 2: The fraction of DIS events with forward photons in the kinematic region 6 < < 
100 GeV^ and 0.05 < y < 0.6 and the pseudorapidity of the photon rj > 7.9. For each 
measurement, the statistical, the total systematic, the uncorrelated systematic uncertainties, and 
the bin-to-bin correlated systematic uncertainties due to the FNC absolute energy scale, the 
impact position of the FNC and the model dependence of data correction are given. 
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Figure 3: The normalised cross sections for the production of forward photons in the pseudo- 
rapidity range r] > 7.9 in deep-inelastic scattering in the kinematic region 6 < < 100 GeV^ 
and 0.05 < y < 0.6 as a function of the longitudinal momentum fraction of the leading 
photon in the range 0.1 < x^l^"^ < 0.7. The data are compared to two predictions of the DJAN- 
GOH Monte Carlo simulation, using LEPTO and CDM to simulate higher orders. Also shown 
are models of hadronic interactions, QGSJET, EPOS and SIBYLL. The lower row shows the 
ratios of the Monte Carlo predictions to the data. The error bars show the total experimental 
uncertainty, defined as the quadratic sum of the statistical and systematic uncertainties. 
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Figure 4: The normalised cross sections for the production of forward photons in the pseudo- 
rapidity range r] > 7.9 in deep-inelastic scattering in the kinematic region 6 < < 100 GeV^ 
and 0.05 < y < 0.6 as a function of the transverse momentum p^"'^ of the leading photon in 
the energy range 0.1 < x^'^"' < 0.7. The data are compared to two predictions of the DJAN- 
GOH Monte Carlo simulation, using LEPTO and CDM to simulate higher orders. Also shown 
are models of hadronic interactions, QGSJET, EPOS and SIBYLL. The lower row shows the 
ratios of the Monte Carlo predictions to the data. The error bars show the total experimental 
uncertainty, defined as the quadratic sum of the statistical and systematic uncertainties. 
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Figure 5: The normalised cross sections for the production of forward photons in the pseudo- 
rapidity range r] > 7.9 in deep-inelastic scattering in the kinematic region 6 < < 100 GeV^ 
and 0.05 < y < 0.6 as a function of the sum of the longitudinal momentum fractions of pho- 
tons xl"™. The data are compared to two predictions of the DJANGOH Monte Carlo simulation, 
using LEPTO and CDM to simulate higher orders. Also shown are models of hadronic inter- 
actions, QGSJET, EPOS and SIBYLL. The lower row shows the ratios of the Monte Carlo 
predictions to the data. The error bars show the total experimental uncertainty, defined as the 
quadratic sum of the statistical and systematic uncertainties. 
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Figure 6: The fraction of DIS events with forward photons in the kinematic region 6 < < 
100 GeV^ and 0.05 < y < 0.6 and the pseudorapidity of the photon 77 > 7.9 as a function of 
and xbj- The inner error bars shows the quadratic sum of the statistical and the uncorre- 
cted systematic uncertainties. Shaded band shows the correlated systematic uncertainties. The 
expectation from the LEPTO and CDM models are also shown. 
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